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Abstract 

The overall characteristics of the solar and atmospheric neutrino oscillation are approximately 
consistent with a tribimaximal form of the mixing matrix U of the lepton sector. Exact tribimax- 
imal mixing leads to #13 = 0. However, recent results from Daya Bay and RENO experiments 
have established a nonzero value for #13. Keeping the leading behavior of U as tribimaximal, 
we use a generalized Fridberg-Lee neutrino mass model along with a complementary ansatz to 
incorporate a nonzero $13 along with CP violation. We generalize this model in two stages: In 
the first stage we assume /r — r symmetry and add imaginary components which leads to nonzero 
phases. In the second stage we add a perturbation with real components which breaks the fi — r 
symmetry and this leads to a nonzero value for #13. The combination of these two generalizations 
leads to CP violation. Using only two of the experimental data, we can fix all of the parameters 
of our model and predict not only values for the other experimental data, which agree well with 
the available data, but also the masses of neutrinos and the CP violating phases and parameters. 

These predictions include the following: (m„ e } « (0.033 — 0.037) eV, (m^) « (0.043 — 0.048) eV, 

(m Vr ) « (0.046 - 0.051) eV, and 59.21° < 5 < 59.34° 

1 Introduction 

The remarkable experimental achievements resulting in experimental data for solar, atmospheric, ac¬ 
celerator and reactor neutrinos pQ give us information about neutrino masses and mixing which can 
be summarized as follows in Table [T] [2J13]: 


Parameter 

The experimental data 

The best fit (±1 <t) 

Am'i^lO^eV 2 ] 

7.12-8.20 

7.43 - 7.81 

\Anrii il \[l§~ i eV 2 ] 

2.13-2.74 

2.21 - 2.64 

2.46-2.61 

2.37-2.50 

sin 2 0 i 2 

0.27-0.37 

0.303- 0.336 

sin 2 0i3 

0.017-0.033 

0.0218- 0.0275 

0.0223- 0.0276 

sin 2 023 

0.36-0.68 

0.37-0.67 

0.400 - 0.461 and 0.573 - 0.635 
0.569- 0.626 

5 

CM 

1 

O 

CM 

1 

O 


Table 1: The experimental data for the neutrinos’ mixing parameters. In the last column, when 
multiple sets of allowed ranges are stated, the upper row corresponds to normal hierarchy and the 
lower raw to inverted hierarchy. 
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The lepton mixing matrix in the standard parametrization is given by mm, 



f C12C13 

S12C13 

Size - 15 \ 

1 1 

f 1 

0 

° \ 

j 

UpMNS = 

^S12C23 — Ci2S23Si3e I<5 

C12C23 — Sl2S23Sl3e l<5 

S23C13 


0 

e ip 

0 

\ 

l S12S23 — Ci2C23Sl3e I<5 

— C12S23 - Si2C2 3 Si 3 e l& 

C23C13 / 

1 1 

^ 0 

0 

e l ° j 


where dj = cos % and s,j = sin (9,;,- (for i, j = (1, 2), (1, 3) and (2, 3)). The phase S is called the Dirac 
phase, analogous to the CKM phase, and the phases p and a are called the Majorana phases and 
are relevant for Majorana neutrinos. However, we should mention that recently the advantages of the 
original symmetrical form of the parameterizations of the lepton mixing matrix has been discussed |5j. 

A relatively successful phenomenological Ansats for the neutrino mixing matrix was proposed by 
Harrison, Perkins and Scott in 2002 [0], and is known as the tribimaximal mixing matrix (TBM). It 
has the following form, 


Utbm 



( 2 ) 


In this matrix the third mixing angle, i.e. 9 13 , is exactly zero. However the results of the Daya 
Bay and RENO collaborations have shown that $13 = 0 is now rejected at a significance level higher 
than 8 cr. A combined analysis of the data coming from T2K, MINOS, Double Cliooz, and Daya Bay 
experiments shows that the best-fit value of 9i 3 is given by sin 2 = 0.026(0.027)1o qo 4 for normal (or 
inverted) mass hierarchy. From these significant data, it has to be concluded that the simple-picture of 
tribimaximal mixing matrix fails. However the smallness of $13 as compared to the other two mixing 
angles encourages us to examine whether a small perturbation on the basic tribimaximal structure 
could lead to a nonzero value for $13 and a more realistic neutrino mixing matrix. 

A successful phenomenological neutrino mass model with flavor symmetry that is suitable for both 
the Dirac and the Majorana neutrinos was proposed by Friedberg and Lee (FL). This model is in¬ 
teresting because when p, — r symmetry is assumed, the resulting neutrino mixing matrix reduces to 
the Utbm- Usually the ansatz for the neutrino mass matrices contain more parameters than can be 
measured in realistic experiments. Thus, the presence of certain conditions or simplifications for the 
neutrino mass matrix is useful. What might first come to mind is the presence of zeros in the mass 
matrix [7]. However, conditions on basis independent quantities, namely the trace and the determi¬ 
nant, are more suitable to consider. The simplest conditions are to set these quantities to zero. The 
condition det(M u ) = 0 [ 8 : leads to zero mass for one of the neutrinos. A zero determinant can be 
motivated on various grounds [9!. The second simplest basis independent condition is a vanishing 
trace, i.e., Tr{M v ) = 0, which is called the complementary ansatz. Its consequences have first been 
investigated in m applying a three neutrino framework that simultaneously explains the anomalies of 
solar and atmospheric neutrino oscillation experiments as well as the LSND experiment. In HU a CP 
conserving traceless M v has been investigated for the simpler case of explaining only the atmospheric 
and solar neutrino deficits. Motivations for traceless mass matrices can be provided by models in which 
Mv is constructed through a commutator of two matrices, as it happens in models of radiative mass 
generation m ■ More interestingly, an approximately traceless M„ can be the consequence of exact 
b — t unification at high energy scales within type II see-saw models m , which in this framework is 
also the reason for maximal atmospheric neutrino mixing Hang. The type II see-saw mechanism was 
the original motivation of the traceless M v condition as investigated in |T0J ■ Also in jT 6 ] the condition 
Tr(M„) = 0 has been used to determine the values of neutrino masses, and the CP phases when CP 
violation is considered. In the present paper we will use the FL model as a starting point, and assume 
the complementary ansatz. 

In the FL model the mass eigenstates of three charged leptons are identified with their flavor 
eigenstates. Therefore neutrino mixing matrix can be simply described by 3 x 3 unitary matrix U which 
transforms the neutrino mass eigenstates to the flavor eigenstates, ( v e , v^,v T ). The Dirac neutrino mass 
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operator in the FL model can be written as 


Mfl = a(v T -v tl ){y T -v ll ) + b{y tl -v e )(v tl ,-v e ) + c(v e -v T ){y e -VT) 

+ mo (v e v e + UftVfj. + v T v T ). (3) 

All the parameters in this model (a, 6, c and Too) are assumed to be real. In the original FL model, also 
known as the pure FL model, mo = 0 and in this case Mfl has the following symmetry v e — > i/ e + z, 
Vp —v Vfj + z, and u T — > v T + z, where z is an element of the Grassman algebra. For constant 2 , this 
symmetry is called the FL symmetry m in which case the kinetic term is also invariant. However 
the other terms of the electroweak Lagrangian do not have such a symmetry. The too term breaks 
this symmetry explicitly. However we may add that the FL symmetry leads to a magic matrix and 
this property is not spoiled by the too term. The magic property has many manifestations which we 
shall discus in detail. Also it has been reasoned that the FL symmetry is the residual symmetry of the 
neutrino mass matrix after the SO(3) x U( 1) flavor symmetry breaking Q2J. The mass matrix can be 
displayed as, 

( 6 + c + too —b —c \ 

—b a + b + mo — a 1 , (4) 

—c —a a + c + too J 

where a oc (Y^ T + Y TI1 ), b oc ( Y eT + Y Te ) and c oc ( Y Te + Y eT ) and Y a p denote the Yukawa coupling 
constants. The proportionality constant is the expectation value of the Higgs field. It is apparent that 
Mfl possesses exact n — r symmetry only when b = c. Setting b = c and using the hermiticity of 
Mfl , a straight forward diagonalization procedure yields U T MflU = Diag m^\ m^} where, 


( 0 ) 

m\ 

= 36 + too 




0 0 

W CN v— 'CO 

= rn 0 , 

= 2 a+ 6 +toq. 

and 

U = Utbm, 

(5) 


Note that in the pure FL model one of the neutrino masses is exactly zero. For a general exact TBM 
neutrino mixing, regardless of the model, the mixing angles are given by sin 2 6*12 = ^ ($12 ~ 35.3°), 
sin 2 023 = \ (023 = 45°), 0 i 3 = 0 and the CP-violating phases can be chosen to be zero. In order to 
have CP-violation in the standard parametrization given in Eq. ©, the necessary condition is S ^ 0 
and 0 i 3 7 ^ 0. In this model these conditions necessarily mandate that jx — r symmetry should be broken. 
Another interesting question is whether 023 = 45° holds after the n — r symmetry breaking. 

There are four independent CP-even quadratic invariants, which can conveniently be chosen as 
^ 11 ^ 11 ,^* 3 ^ 13 ,^ 21^21 and U 23 U 20 and three independent CP-odd quartic invariants [19], 

J = ImiUuU^M 

h = lm[(U^U 12 ) 2 ] 

I 2 = lm[(U^U 13 ) 2 ). ( 6 ) 

The Jarlskog rephasing invariant parameter J EDI, is relevant for CP violation in lepton number 
conserving processes like neutrino oscillations, ii and I 2 are relevant for CP violation in lepton number 
violating processes like neutrinoless double beta decay. Oscillation experiments cannot distinguish 
between the Dirac and Majorana neutrinos. The detection of neutrinoless double beta decay would 
provide direct evidence of lepton number non-conservation and the Majorana nature of neutrinos. 
Many theoretical and phenomenological works have discussed massive neutrino models that break 
p — r symmetry as a prelude to CP violation[211. 

In this paper we start with the simple FL model with exact fi — t symmetry, which leads to Utbm 
mixing matrix. We then generalize this model by adding complex parameters to the elements of Mfl 
which can ultimately be linked to complex Yukawa coupling constants. This generalization can be 
broken down to two clearly distinguishable pieces. In this paper we first study the results of each piece 
separately and then investigate the results when both pieces are applied simultaneously. Both pieces 
are separately exact solvable, while the combination is not. First we add a real term proportional to the 
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matrix representation of the permutation group element S 132 , merely to break the /j, — t symmetry. As 
we shall show, it suffices to treat this piece as a perturbation. We find that this generates a nonzero 0 13 
whose smallness justifies the use of perturbation theory. Second we add arbitrary imaginary coefficients 
while preserving /x — r symmetry. This results in the generation of nonzero phases, including 5 , which 
turn out to be large. For combining the two steps, we first do the complexification step exactly and 
finally add the fi — r symmetry breaking perturbation. We then find that nonzero values are generated 
for both 0 i 3 and <5, and this leads to CP violation. 

This paper is organized as follows. In section 2, we introduce our model and show the results 
of the two aforementioned generalization, separately. We then combine the two generalization and 
show the results. We use some self consistency arguments along with the complementary ansatz to 
reduce the number of our free parameters. In section 3, we map two of the experimental data onto 
the allowed region of our parameter space. We then check the consistency of all other experimental 
data with the two overlap regions, which selects only one region. This region is extremely small and 
almost pinpoints all of our parameters. The results indicate that the complexification had to be treated 
nonperturbatively, since the phases turn out to be large. However the coefficient for fj, — r symmetry 
breaking justifies the use of perturbation theory for that part, since $13 turns out to be small. We 
then not only check the consistency of all of our results with the available experimental data, but also 
present our predictions for the actual masses and CP violation parameters. In section 4, we state our 
conclusions. 


2 The Model 

In this section, we study the effects of generalizations of the FL model that break the /r — r and CP 
symmetries. In the first stage we perturb the mass matrix so as to break the /i — r symmetry. We 
choose the mass matrix to be of the form M v = M 3 + M„ = M® + \S Z , where M 3 = Mfl with b = c, A 
is a perturbation parameter with dimension of mass, and S l is an element of the permutation matrices 
which breaks the fi — r symmetry. For example for S l = S^ 132 \ we have, 


M v 


M° + Af' = M° + AS (132) 


/ 2 b + mo 

-b 

-b \ 

/ ° 

1 

0 \ 

-b 

a + b + mg 


+ A 0 

0 

1 

\ —b —a 

a + b + too / 

' V 1 

0 

0 / 


(7) 


Notice that M v and Mfl are both magic and symmetric matrices since they both commute with the 
magic S matrix defined by 


S = 


F 

T 

T 

T 

F 

T 

T 

T 

F 


( 8 ) 


Therefore one of the eigenstates must be (-^, -^). The elements of the basis in which is 

diagonal are simply the columns of Utbm , as shown in Eq. ©, and are as follows, 



(9) 


Here M 3 and M v are real matrices. Needless to say, this may generate a non-zero $13 but will not 
lead to CP violation, since this computation necessarily yields <5 = 0. As stated in the Introduction, 
although this problem is exactly solvable it not only suffices to treat this part as a perturbation since 
013 turns out to be small, but also it is advantageous to do so for comparison purposes with the part 
when all generalizations are combined. The perturbation expansion of the mass eigenstates of M v to 
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first order is 


,( 0 ) 


*7U 


(v 0 ) \m'M 0) ) 


+ ~ v j 


(o) 


The Vi& are simply the columns of the mixing matrix U and are as follows, 


( 10 ) 


v\ 


( 


\ 



(a-b) 




( V3 ^ 


( 2+2 (a-b) \ 

^2 = 


) ^3 = 

1 4. * 

+2 ' 4+2(a—b) 


\ V3 / 


\ s/2 + 4s+2(a-b) ) 


( 11 ) 


Comparing 1/3 given in Eq. ED with the third column of Upmns as given by Eq. m, we immediately 
obtain all of the mixing angles ($ 13 , $ 12 , $ 23 ) and the CP violating phase 6 in terms of the parameters 
of our model, i.e. a , 6 and A, as follows, 


sin 2 0 i 3 


sin 2 0 i 2 


sin 2 023 
tan 5 


2A 2 

16(a — b) 2 + 3A 2 ’ 

1 16(a — b) 2 + 3A 2 

3 cos 2 0 i 3 48(a — b) 2 + 3A 2 ’ 

1 4A(a — b) 

2 ~ 16(a - b) 2 + A 2 ’ 

0 . 


( 12 ) 


Notice that we managed to obtain a nonzero 0 i 3 . 

The first order corrections to the neutrino masses are obtained from = m^Sij. 

Therefore the masses up to first order correction are as follows, 

=3 b + toq — —, m ^ = mo + A, m 3 ^ = 2 a + b + mo — —. (13) 


From the complementary ansatz we obtain, a = ~( 4 b + 3 m o) _ j n limit A —^ 0, the mixing matrix U 
must reduce to Utbm- Therefore we find (a — b) < 0, which has also been obtained in [ 22 | . Using 
these conditions, along with the usual convention of FL model in which mo > 0, one can easily show 
that — ^mo < 6 < 0, — |mo < a < —^?no, m 3 < 0, m 3 < mi and Am§i = m 3 — m\ > 0. Therefore, 
we have normal hierarchy at this stage. As we shall see this property remains true in all stages of our 
model. 

In the second stage we let all of the coefficients a, 6 , c of the original Mfl be complex with the 
restriction b = c which preserves the exact n — r symmetry of the original M°. This mass matrix can 
be displayed as 


/ 2(b + iB)+mo —(b + iB ) —(b + iB ) \ 

M£° = I —(b + iB) (a + iA) + (b + iB)+mo —( a + iA ) I, (14) 

y — (b + iB) —(a + iA) (a + iA) + (b + iB) + mo ) 

Where the superscript “c” stands for “complex”. The M£° is a non-Hermitian matrix, so in general 
we need two distinct unitary matrices U and V to diagonalize it, i.e. M^ ag = WM^°V. These 

matrices can be easily obtained by diagonalizing and M£°, separately. The matrices U 

and V are the conventional transformation matrices for the left-handed and right-handed neutrinos, 
respectively. An interesting point is that for one can easily show that U = V = Utbm, due to 
the /I — t symmetry. Therefore the basis in which M“° is diagonal is the same as that of M®, i.e. 
{Vi°\ i = 1,2,3} given by Eq. ([21). The eigenvalues of M“° are as follows, 

mf = (36 + mg) + i(3B), mf = mo, = (2a + b + mo) + i(2A + B). (15) 
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From the complementary ansatz, i.e. Tr(M£°) = 0, we obtain, 


-(46 + 3?7i 0 ) 


and A = —2 B. 


( 16 ) 


It is worth mentioning that in this case because of the imaginary terms in the mass matrix we 

have two phases that appear in the mass eigenvalues shown in Eq. (USD , i-e. exp[iarctan( 35 ^_)] and 
exp[?’ arctan( 2 a+^n 0 ^' ^or Dirac neutrinos these phases can be removed and for the Majorana 
neutrinos these phases remain as Majorana phases and contribute to CP violation. 

Notice that for the complex FL model with /j, — r symmetry 6 13 = 0. Moreover, the “zero sum” 
condition, i.e. mf 1 + m 2 ° + mf = 0, is equivalent to the traceless condition, since Tr(Mff) = 
Tr(WVMfjff ag ) = Tr(M^ ag ) = 0 . However, notice that in general Tr(M) ^ Tr{M < n a g) [II] - 

The experimental data have now definitely confirmed that Ajtj^ = to| — m\ > 0. Therefore, using 
Eq. (fT51) , we can define the following new variable, 


B 


\y\ 



(17) 


where — 1 < y < 1 and we have defined dimensionless parameters: 6 = — B = — and a = —. 
The limit \y\ —> 1 leads to |to 2 °| —> The conditions stated for a and 6 below Eo. (fTdl) . i.e. 

— i < 6 <0 and — | < a < — |, are consistent with the definition of B . 

In the final stage we combine two basic generalization to the Mfl, namely the addition of imaginary 
components while preserving y — r symmetry as a finite distortion, and the addition of a permutation 
matrix with a real coefficient to break the y — r symmetry as a perturbation. Then we obtain as 
follows, 


M?. 


m: 


cO 


M„ 


2(6 + iB) + too 
-(6 + iB) 

— (6 + iB) 


-(6 + iB) 

(a + iA) + (b + iB) + m 0 
— (a + iA) 


-(6 + iB) 

— (a + iA) 

(a + iA) + (6 + iB) + m 0 


/ 0 1 0 
A 0 0 1 
\ 1 0 0 


(18) 


where is a complex symmetric non-Hermitian matrix given in Eq. (1141) and M' = A S<- 132 \ as 
given by Eq. 0, is real matrix. The columns of the mixing matrix U in Eq. 0 are eigenstates of 
III/ M£ 0 + M£°*M V + M'JM£°, where we have dropped a term which is 0( A 2 ). To proceed, 
we recall that Mff is Hermitian and its eigenstates are the same as the columns that produce 
Utbm, and its eigenvalues are |m^ c0 ^| 2 , Im^) 2 , and Im^ 0 ^ 2 . The basis in which is diagonal is 
given by, 


,(°) 


E 


{ mf *) 2 — ( m ^ 0 ) 2 


,(o) 


(19) 


6 







A straight forward calculation yields, 


v 


c 

1 


2 v / 3 



(36B 2 + (3-2m 0 ) 2 )A 2 
(l+2h) 2 (3-2m 0 ) 2 



(3—2mo — i6B)\ 
2v^(l+2b)(-3+2m 0 ) 

(3—2mo—i6B)\ 


\ 

(?) 

■v'i = 

73 

) 

V 73 / 


( 


\ 


— \Z^"(3—2mo+^6B)A 

(l+2b) (-3+2m 0 ) 

1 (3-2m 0 +i6B)A 

v^ i ~6y/2(l + 2b)(-3+2m n ) 

1 /, , (36B2 + (3-2m 0 )2)A2~ 

2 V H (l + 26) 2 (3-2m 0 ) 2 

)_ (3-2m n +.6B)A 

75 6v2il + 26l< 3+2m„ 

1 /, , (36B2 + (3-2m 0 )2)A2~ 

2 V H (1 + 26)2 (3 —2mg)2 


\ 


/ 


( 20 ) 


These are also simply the columns of mixing matrix, U, that contain CP violating phases. Comparing 
Eq. (0(71) with Eq. dTJ), we immediately obtain all of the mixing angles (0i3,0i2,023) and the CP-violating 
phase in terms of 6 , B and A = — as follows, 


2(1 + 36S 2 )A 2 

36(1 + 2b') 2 + 3(1 + 36H' 2 )A' 2 ’ 

24(1+ 26') 2 

1 ~~ 36(1 + 2b') 2 + (1 + 36B' 2 )A' 2 ’ 

1 , 6(1 + 2 b')\' 

2 + 36(1 + 2b 1 ) 2 + (1 + 36fT 2 )A' 2 ’ 

6 B . (21) 

The hrst order corrections to the neutrino masses are obtained from M v + 

M^Mf\vf ] ). Therefore by using Eq. (fT5l) and the first order of neutrino mass correction we have, 

= |m ° 0 | 2 - (36 + m 0 )A, 

= |m. 2°| 2 + 2 m 0 A, 

= | m.|° | 2 — (2a + 6 + mo)X. (22) 

As in Eq. (1151) . and m 3 are complex. Extracting the real and imaginary components of the masses 

and introducing them as two phases in accordance with Eq. CD, we obtain p = — arctan ( 3b+ ^ _a ) an d 

a = — arctan ]j + arctan ( 2 a+bVm 0 - A ) ■ Since M£° is a symmetric matrix, it could also be 
used as a Majorana mass matrix, and these phases can be considered as the Majorana phases. In 
the Dirac case these phases are transferred to the mass eigenstates. However in the Majorana case 
the phase factors remain and contribute to the CP violation. The requirement det(U) = 1 leads to 
p + a = 0 [23 . This in turn implies the following, 

A' = 2 (V + y / '(1 + 26') 2 + 3H' 2 ^ . (23) 

Using Eq. (1221) , we can write the neutrino mass squared differences as follows, 

A rri 2 i = |m 2°| 2 — |m 3 °| 2 + 3 (6 + m 0 ) A, 

Am 2 i = K°| 2 -K°| 2 -2(a-6)A, 

Ato 32 = |mf| 2 -|m^°| 2 -( 2 a + 6 + 3m 0 )A. (24) 

A rephasing-invariant measure of CP violation in neutrino oscillation is the universal parameter J 
[2D] given in Eq. JD]), and it has a form which is independent of the choice of the Dirac or Majorana 


sin 2 0 i 3 = 

sin 2 0 i 2 = 

sin 2 023 = 

tan S 
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neutrinos. Using Eq. EU1) the expression for J simplifies to, 

4(l+ 26') B A' 

J “ 12 + 486' (1 + b ') + (1 + 36B' 2 ) A ' 2 ' 

Notice that we have reduced the number of free parameters of our models to just two, i.e 
using all the self consistency arguments and the complementary ansatz. 

3 Comparison with experimental data 

In this section we compare the experimental data with the results obtained from our model. We do 
this by mapping two of the constraints obtained from the experimental data onto our parameter space, 
b' and y : as shown in Figure ©■ The two most restricting experimental data come from the values of 

and sin 2 $ 13 . The overlap of and sin 2 $13 and our model is restricted to two tiny regions 

close to the top left corner of the parameter space. 


(25) 
. b and y, 



b 


Figure 1: In this figure the whole region of the \y\-b plane which is allowed by our model is shown. The 
overlap region of the experimental values for and sin 2 #13 with our model are two tiny regions close to 

Am 31 

the top left corner of the rectangular. In the zoomed box we have magnified the overlap regions. As mentioned 
in the text, only the lower overlap region is consistent with all of the experimental data. 


First we check the consistency of the rest of the experimental data with the overlap region. We 
find that the upper overlap region is barely ruled out by the experimental value for sin 2 623 , and the 







































lower overlap region is consistent with all the rest. Having determined the unique overlap region, we 
can not only determine all of the parameters of our model, but also predict the masses of the neutrinos 
and the phases. The results for the parameters are as follows, 


mo « (3.3 — 3.7)10~ 2 eP, 
a « -(1.71-1.93)10 _2 eP, 
b « -(1.61-1.81)10“ 2 eP, 

B w (0.92 — 1.04)10 _2 eP, 

A w -(1.98- 2.96)10“ 4 eP (26) 

By comparing the values of these parameters we conclude that the complexification part of the general¬ 
ization of Mfl , accomplished by adding imaginary components denoted by B, could not be treated as 
a perturbation, while the second part of the generalization, accomplished by adding AS 1123 , can indeed 
be treated as a perturbation. This is precisely what we have done. 


Parameter 

The exp. data 

The best exp. fit(±lcr) 

Predictions of our model 

Am^lO^eP 2 ) 

(7.12-8.20) 

(7.43-7.81) 

(7.12-7.84) 

ATO 3 1 (10 3 eH 2 ) 

(2.26-2.77) 

(2.46-2.61) 

(2.23- 2.80) 

sin 2 6*12 

0.27-0.37 

0.30- 0.336 

0.331- 0.333 

sin 2 0 i 3 

0.016-0.037 

0.0281 - 0.0275 

0.0218- 0.0275 

sin 2 023 

0.39-0.64 

0.400 - 0.461 

0.573 - 0.635 

0.439- 0.446 

S 



59.21° - 59.34° 

a = -p 



59.44° - 59.51° 

\A 



0.029- 0.033 

masses 



mi « (0.031625 - 0.035687)eH, 
(m Ve ) « (0.033- 0.037) eH 



m 2 « (0.032732 - 0.036769)eP, 
(m„ u ) w (0.043-0.048) eP 



m 3 w (0.056849 - 0.063819)eP, 
{m„ T ) « (0.046- 0.051) eP 


Table 2: The available experimental data for neutrinos for the case of normal mass hierarchy and the 
predictions of our model. These predictions are obtained from our parameters as shown in Eq. (1261) . 
Specifically we have used b « —(0.489 — 0.4904) and y ~ (0.948932 — 0.956141). 


Having determined all of the parameters of our model, we can compare our results with the exper¬ 
imental data. In Table [2] we state all of the relevant experimental data presented at 3<r [2] , along with 
the predictions of our model. As shown in Table [5] we have predictions for some physical quantities for 
which no experimental data exist. We should mention at this point that if we had chosen S ' 123 instead 
of S 132 to break the fj, — r symmetry and generate a nonzero value for # 13 , all of our final results and 
prediction would be unaltered except for J —> — J. 

An important experimental result for the sum of the three light neutrino masses has just been 
reported by the Planck measurements of the cosmic microwave background (CMB)[24j, which is 

Y rn v < 0.23eF(P/ancfc+WP+highL+BAO), (27) 

This sum in our model is ]T) m v ss (0.12 — 0.14)eP, which is consistent with the above constraint. For 
the flavor eigenstates only the expectation values of the masses can be calculated and they are obtained 
by the following relation, 

3 

(m Vi ) = Y \ Ui o\ 2 \ m jl ( 28 ) 

3 =1 
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where i = e, p, and r. Our predictions for these quantities are also shown in Table [2] The Majorana 
neutrinos can violate lepton number, for example in neutrinoless double beta decay (0^/3/3) [25J. Such a 
process has not yet been observed and an upper bound has been set for the relevant quantity, i.e. (m J/ee ). 
Results from the first phase of the KamLAND-Zen experiment sets the following constraint (m Vee ) < 
(0.12 — 0.25) eV at 90% CL [26]. Our prediction for this quantity is {rn„ ee ) ~ (0.032 — 0.036) eV which 
is consistent with the result of kamLAND-Zen experiment. We also predict ( m Ve ) ~ (0.033 — 0.037) eV, 
(m v J w (0.043 - 0.048) eV and (m„J « (0.046 - 0.051) eV. 

4 Conclusions 

In this paper we propose a generalized Friedberg-Lee neutrino mass model in which CP violation 
is possible. Our generalization consists of two steps. First we add imaginary components to the 
mass matrix, such that the resulting matrix is symmetric, magic with zero trace. The results is that 
p — r symmetry is preserved and nonzero phases including S are obtained. However, sin $13 remains 
zero. In the second stage we break the p — r symmetry by adding a real matrix which is the matrix 
representation of the permutation group element S' 1 - 123 ). This part is treated as a perturbation and 
a nonzero $13 is generated. The combination of these two generalization then produces CP violation. 
We use self consistency arguments and the complementary ansatz to reduce our parameter space to 
two dimensions. Mapping two of the experimental data, i.e. the allowed ranges of sin 2 $13 and , 
onto our allowed region of the parameter space almost pinpoints the values of our parameters. We 
then check the consistency of all other experimental data with the allowed range. We then predict the 
mixing angles, the masses of neutrinos, the CP violation parameters <5, p, <7, and J. All of our findings 
are shown in Table [2] along with the latest experimental values. As is evident from Table [2j all of 
our predictions are within the allowed ranges determined by the experiments and there is an excellent 
agreement between these values. Our predictions for the neutrino masses and CP violation parameters 

5 and J are to tested by future experiments. 
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